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Abstract 

The antiferromagnetic Heisenberg model on a chain with nearest and next 

nearest neighbor couplings is mapped onto the SO(3) nonlinear sigma model 

in the continuum limit. In one spatial dimension this model is always in its 

disordered phase and a gap opens to excited states. The latter form a doubly 

degenerate spin-1 branch at all orders in 1/N. We argue that this feature 

should be present in the spin-1 Heisenberg model itself. Exact diagonaliza- 

tions are used to support this claim. The inapplicability of this model to 

half-integer spin chains is discussed. 
75.10.Jm, ll.10.Lm 
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I. INTRODUCTION 



Interest in low-dimensional magnetic systems has been great in recent years, partly be- 
cause of the widespread belief that magnetism plays a key role in high-temperature supercon- 
ductivity, but also because of the successful application of field-theoretic methods to these 
systems, in particular to spin chains. Indeed, a mapping from the spin-s antiferromagnetic 
(AF) Heisenberg chain with nearest-neighbor (NN) coupling to the 0(3)/ 0(2) sigma model 
has led Haldane to conjecture^ that integer spin chains should exhibit a gap to a triplet 
of excited states, whereas half-integer spin chains should not. This conjecture has later 
been confirmed by numerical calculations JlH Such a mapping may be termed semi- classical, 
since it is constructed by introducing a local field characterizing the order expected in the 
classical ground state and by taking into account fluctuations of the spin variables around 
this local order. In the case of the simple AF chain, a collinear order exists in the classical 
ground state and the local field introduced is a unit vector e(x, t) representing the staggered 
magnetization. The long wavelength effective action obtained for this field is that of the 
0(3)/ 0(2) sigma model, with Lagrangian density 



The coupling g and the velocity v are related to the spin s, the AF coupling J and the 
lattice spacing: v = 2 J as and g = 2/s. To this action one must add a topological term S top 
for half-integer spin (cf. Sec. IV). An extension of this mapping to dimensions higher than 
one has been obtained by many authors, with the difference that no topological term arises, 
hence no distinction between half-integral and integral spin. 

D ombre and ReaJ have conducted a similar analysis for the antiferromagnetic Heisen- 
berg model on a triangular lattice. The essential difference here is the presence of frustration, 
leading to a classical ground state characterized by a 120° order. The local order must then 
be specified by a rotation matrix instead of a unit vector and the order parameter is thus 
an element of SO (3). The long- wavelength action found in this case is the 5*0(3) nonlinear 
sigma model, with Lagrangian density 




(1) 
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C=—{~ix(d t R-%R)-vtT(PVIC 1 -VR)} . (2) 

Here -R(x, t) is a position- and time-dependent rotation matrix and P is the constant diagonal 
matrix diag[l, 1,0]. Again, the constants g and v depend on the spin s, the lattice spacing 
a and the AF couplings J. At zero temperature and in dimension two, this model has an 
ordered phase for g < g c and a disordered phase otherwise. Since g ~ a/s, the spin-| 
case is the closest to the disordered phase. However, it is now widely believed! that the 
ground state of the spin-| antiferro magnet on a triangular lattice has long range order and 
consequently the disordered phase of the SO (3) model is not physically realized to two- 
dimensional antiferromagnets, at least at half- filling (one spin per site). 

In this work we will argue that this disordered phase of the SO (3) sigma model could be 
realized in a frustrated antiferromagnetic chain, with Hamiltonian 

H = J ^ S n • S n+ i + J' ^ S n • S n+ 2 • (3) 

n n 

Introducing a next-nearest-neighbor (NNN) coupling J' modifies the classical ground state 
if J' / J > 4: the spins still lie on a common plane, but instead of being collinear (as in 
the Neel state), they are arranged in a canted configuration (see Fig. |1|) in which each spin 
makes an angle a with its predecessor, given by 

cosa = ~~U> ' ^ 
In the special case J' / J = |, the classical order is quite similar to that of the triangular 
lattice, with its 120° angle from site to site and its periodicity of three sites. For spin-| 
this particular case constitutes the Majumdar-Gosh model, whose ground state is exactly 
knownE For a generic value of J' / J the classical order is incommensurable, with infinite 
periodicity. The spin-1 case has been studied previously using various methods: exact 
diagonalizations@, approximate mappings to field theories involving fermionsil or bosonJ^, 
among others. 

This paper is organized as follows. In Sec. II we explain how to obtain the 5*0(3) 
nonlinear sigma model as the continuum limit of the model defined in Eq. (H). In Sec. Ill 



we argue that the disordered phase of this model is characterized by a singlet ground state 
with a gap to two degenerate triplets of excited states. Hence this model cannot adequately 
represent the half-integer NNN chain, but may be correct for integer spin. In Sec. IV we 
discuss our results with the help of exact diagonalization results and discuss the apparent 
lack of distinction between integer and half-integer spin in this model. 

II. SEMI-CLASSICAL MAPPING TO THE SO(3) NONLINEAR SIGMA MODEL 

In this section we argue that the long wavelength, continuum theory describing the 
frustrated antiferromagnetic chain is specified by the Lagrangian density of Eq. ([17|), with 
the notation explained thereafter. 

The Lagrangian description of spin dynamics, as used in path integrals, requires the 
introduction of spin coherent states. The unfamiliar reader is referred to Fradkin's text0 or 
to Manousakis' review0. Each quantum spin Sj is described by a fluctuating unit vector rij 
and the action associated to a spin Hamiltonian H is 



wherein Sj is replaced by srij in the Hamiltonian. A is the vector potential of a magnetic 
monopole of flux 4tt. Thus, given a closed curve n(£) on the unit sphere, the integral J A - dn 
equals the area of the sphere enclosed by the curve, modulo 47r. 

We shall use a description of each spin Sj in terms of a slowly varying 5*0(3) order 
parameter R(x,t) and a local magnetization vector l(x,t): 



Here n£ is the orientation of the spin in a classical ground state taken as reference, in 
which all spins lie on the plane defined by two mutually orthogonal unit vectors i and j; a 
is the pitch of this ground state, as given in Eq. finally, a is the lattice spacing: It is 




(5) 




(6a) 




(6b) 
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assumed that the spin configurations that contribute significantly to the path integral are 
locally close to the classical ground state, and the approximation that the deviation al is 
small will be controlled at the same time as the continuum approximation. 

Some comments are in order concerning this representation. If J' = J/2 (or a = 2n/3), 
the periodicity is 3 and the relation ([]) may be considered as a bona fide change of variables 
if we group the spins in sets of three and assume that the fields R and 1 do not vary within 
such sets. We check that the number of degrees of freedom match: 6 per set of three spins. 
However, it is more convenient to assume that R and 1 vary slowly from site to site, in which 
case the representation (||) may not be regarded as a change of variables, but simply as long 
wavelength description of the fluctuations around the classical ground state; then we need 
not restrict ourselves to the case J' = J/2. 

The next step is to substitute the representation (||) into the action fl5|) and to Taylor 
expand when needed in order to get a continuum action in terms of R and 1 only. We will 
treat the kinetic term first, and then the Hamiltonian. For the sake of convenience, we will 
assume a periodicity of iV spins in the classical ground state. The case N = 3 (or J' = J/2) 
is the simplest, leading to a calculation almost identical to that of Ref. [|. The next simplest 
is iV = 5, with J' « 0.809J, and so on. 

An explicit calculation of the kinetic term is not needed here, since it would be almost 
identical to that of Ref. |], to which the reader is referred. The difference here lies in the 
number of reference vectors n?, which is N instead of 3. The kinetic part of the Lagrangian 
density is then 

C K = as(Tl) • V (7) 

with the definitions 
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V a = -~eabc(R-%R)i>c (8b) 
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wherein the indices in spin space correspond to the axes defined by T, j and i x J, and (• • •) 
means an average over N contiguous sites. 

In order to write down the continuum limit of the Hamiltonian H, we need to calculate 
the interaction of a spin Sfc with its four neighbors and then average the result over 
the N spins S k of the period. Specifically, the Hamiltonian density is 



n 



1-5>M = !<#(*>> . (9) 



2aN Y 2a 



The contribution may be written as follows: 

( 3' 3' } 

Jj{k) _ j J g fc . Q k+l _|_ S fc . Sk-l + ~~jSk ■ Sk+2 + -J Sfc • Sfe_ 2 f 

[ 23' 23' 1 

+ aJ < Sfc • C^Sfc+i — Sfc • C^Sfc-l H ' J~^ k ' ^Sfc-2 [ 

+ ^ |s fc • cPS fc+1 + S fc • «9 2 S fc -i + ^S k ■ d 2 x S k+2 + ^S fc • «9 2 S fc _ 2 j + • • • (10) 

where the spin S n stands in fact for its expression in terms of R and 1 (Eq. ([])) and where 
the derivatives act on these same fields. We will carry the expansion to second order only, 
enough to yield a non trivial Lagrangian. We shall also use a Taylor expanded version of 
the representation (||): 

S fc = sR(n c k ) + saR{\ - [1 ■ n c k ]n c k ) + sa 2 R(n c k ) [§(n£ • l) 2 - ±1 2 ] - sa 2 R(l)[n c k ■ 1] + C(a 3 ) . 

(11) 

We may then write the Hamiltonian density as an expansion in powers of a: 

n = n {0) + an [1) + a 2 n {2) + ■■■ (12) 

7v-°> is a constant, independent of R and 1. Ti.^ vanishes for the following reasons: Part of 
it depends on 1, but it is proportional to (n|), which is zero. Another part depends on R, 
and is given by 

n(1) = ^flM^W - KX-i, 6 » + 2 -AKX +2)6 > - KX_2, 6 })} • 

(13) 



But one easily calculates that 

\ n k,a n k+n,b) 



( cos na sin na \ 



sinna cosna . (14) 
V OJ 

After substituting and using the relation (f|), we find out that the expression ( |i~3|) vanishes. 

Finally, we are left with the second order part 7i^ 2 \ It is a straightforward exercise to 
show that 



U = J -^-hM ab l b - ^ - ±) tr[Pd x R^d x R] 
where £ = — cos a = J / 4 J' and where the diagonal matrix M has elements 

(i + 2 



(15) 



M: 



t 



M ll = M 22 = M^{e-i + \) . 



(16) 



Again, the matrix P is the projector onto the plane defined by T and j, as in Eq. (Q). The 
full Lagrangian density is then C = Ck — 7i. 

The last step consists in integrating out the field 1. Since the latter appears quadratically 
in £, this operation amounts to substituting into C the solution of the classical equations of 
motion for 1. The final result for the Lagrangian density is 



C = Y ^-MQd t R- l d t R] - ctr[Pd x R- l d x R]} 



(17) 



where the matrix Q is diagonal, with elements 

Qu = Q22 = 1 Q33 
and where the constants g and c are defined as 



2^(1-0 

e+(t- 1) 2 



(18) 
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2 1 + ^ 



c = Jas^—^- yj\ - £ 2 . 



(19a) 
(19b) 



In the special case £ = ~, which corresponds to J' = J/2 and a = 27r/3, one recovers precisely 
the form of the Lagrangian @ since Q33 = 1. The parameter £ ranges from (J 7 — ► 00) to 



1 (J' = J/4). The characteristic speed c diverges as £ — > 0, whereas the coupling constant g 
diverges as £ — > 1, the collinear phase boundary. In all cases, the Lagrangian (|T7|) does not 
possess Lorentz invariance, which makes it a qualitatively distinct theory from the versions 
of the 50(3) nonlinear sigma model studied in the context of classical critical phenomen 
Its symmetries comprise global left rotations R — > UlR with Ul G 50(3) and global right 
rotations R — > RUr with Ur G 50(2) (i.e. Ur commuting with P and Q). The first of these 
reflects invariance under rotations in spin space, whereas the other will have consequences 
on the spectrum of the theory: a degeneracy of the excitation branches. 

It is important to stress here that in going from the discrete Heisenberg Hamiltonian @ 
to the continuum Lagrangian fll7|) we have assumed that smooth configurations dominate the 
path integral. There seems to be no distinction between integer and half-integer spins within 
this mapping. In particular, there is no topological term of the type arising in the collinear 
antiferromagnetic chain (J' < J/4). In fact, none could exist, since the relevant homotopy 
group is trivial: 7r 2 (50 (3)) = 0. However, there are qualitative differences between integer 
and half-integer spin in this system. The question as to why they do not appear explicitly 
in the field theory will be discussed in Sec. IV. 



III. SPECTRUM OF THE NONLINEAR SIGMA MODEL 

In this section we discuss the spectrum of the model defined in Eq. ([T7D, in particular 
regarding the gap to excited states and the degeneracy of the excitation branches. To this 
end it is preferable to rewrite the Lagrangian (|17]) in terms of two orthonormal vectors ei 
and e 2 specifying the rotation matrix R. Defining e 3 = e x x e 2 , the elements of the matrix 
R may be expressed as R a b = (e b ) a and trlQdtR^dtR] = J2 a Qaa(d t e a ) 2 . Using the fact that 



(d t e 3 ) 2 = (<9 tei ) 2 + {d t e 2 f - 2( ei • d t e 2 ) 



(20) 



we may finally express the Lagrangian in the following form: 



1 fl 



(d tei ) 2 + (d t e 2 ) 



^(e 1 -^e 2 ) 2 -c[(9 :c e 1 ) 2 + (^e 2 ) 
c L 



(21) 



wherein the coupling constants g, 7 and the velocity c are now 



_ 2 (1 + 0(2^-2^ + 1) 

f ^ (22a) 

1 + e 



c = /as— ^C 1 - e)(2e - 2£ + 1) (22b) 
7 = 4£(1 - (22c) 

The functional integration measure associated with the fields ei j2 must incorporate the 
constraints e\ = e| = 1, e x • e 2 = 0. This may be done by introducing three Lagrange 
parameters an, a 12 and 012 = 021 > and by adding the following constraint term to the 
Lagrangian: 

^■constr = 2^ a b(.^a ' &b $ab) ■ (23) 

Since this model is defined in one spatial dimension, the Mermin- Wagner- Coleman the- 
orem applies and the global SO (3) symmetry is not broken. If we extended the model to 
dimensions greater than one, a broken symmetry phase could be realized, for g below some 
critical value. In such a phase, three Goldstone modes would appear, corresponding to the 
three parameters of the broken SO (3) symmetry. These three Goldstone modes would have 



a linear dispersion relation, two of them with speed c and a third with speed c/Jl — 7/2. 
The most elegant way to see this is to go back to the form flTTD of the Lagrangian and to 
substitute, in the small oscillation approximation, R ab « 6 ab + e abc VL C) keeping the terms of 
order Vl 2 c only. Since the field fl c is unconstrained, the dispersion relations are easily read 
off. In the case of an antiferromagnet on a triangular or hexagonal lattice (£ = |), the two 
speeds are c and ca/2. This may also be observed within spin-wave theory, since it is a 
feature of the ordered phase. 

However, we are interested here in the disordered phase of the model, which we shall 
study within the large- iV approach (here iV is the number of components of the vectors e a , 
normally 3). In the imaginary-time formalism, the partition function of the model is 

Z = J [dei] [de 2 ] [da ab ] exp — / dxdr C E (24) 



with the Euclidian Lagrangian density 
Ce = o { [(^ e i) 2 + (^ e 2) 2 ] - <?7(ei • d T e 2 f + [(^ ei ) 2 + (d x e 2 ) 2 ] } + a ah (e a ■ e b - -5 ab ) 
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(25) 



where we have rescaled e a by a factor yfg in order to recover the standard normalization 
for the kinetic term. The characteristic speed c has been set to unity in order to lighten 
the notation; it may be restored by dimensional analysis. We will also use a Hubbard- 
Stratonovich decomposition of the quartic term: 



cxp 



dxdr ^{e l -d T e 2 f 
2 



J [d(p} exp — J dxdr ||0 2 — v /^70(ei • d T e 2 



(26) 



We then proceed to find the large- iV saddle point. In other words, we assume that 
the auxiliary fields and a ab take a constant value, which is determined by extremizing 
the effective potential obtained by integrating the fields e a . This exercise is better done in 
Fourier space, in which the Euclidian action with constant auxiliary fields may be written 



as 



du dk 
2tt2tt 



e* a (u,k)K ab e b (u,k) + ±L 2 



&ab&ab 
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(27a) 



K(u,k) = 



cu 2 + k 2 + (Tn o- 12 - iuKpy/g^'* 
_ a 12 - iuj(f)^/(py J 1 + k 2 + a 22 



The effective potential is then 

i ,9 1 / \ N f du dk , , T w , s 

Veff = \<P 2 ~ + <r 22 ) + - J — — logdetAT(u;,fc) 

where we assume the system to be limited by a box of side L. 
In terms of the variables 



° = §(cn + 0-22) and rf = \(a u - a 22 ) 2 + a 



(27b) 



(28) 
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(29) 



the saddle-point equations are the following: 
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dV *ff _ n _ 1 N r dujdk 2(uj 2 + k 2 + a) 

u - ~Z + TT / /. o ; ,9 ; _\o — 15 — 17715" .9 ( dUa 



9(7 5 2 y 2vr 2vr (cj 2 + A; 2 + cr) 2 - if - c/7(/> 2 a; 

<9V e / f n f dto dk 1 



<9?7 N V J ^ ^ ^ 2 + fc2 + ^ 2 _ ^ 2 _ ^2^2 ( 30b ) 

dV e ff f duj dk ijj 2 . , 

— = = 0- Ncpgry J -- - — - — — (30c) 

The immediate solution to these equations is = 77 = and a 7^ 0, determined by the 
simpler equation 

I _ N r duo dk 1 

~g~ ^2 J 2^2^ cj 2 + P + ct ^ ^ 

It can be shown explicitly^ that the solution <fi = r\ = is the only acceptable one in this 
case. Eq. (|3lD coincides with the saddle-point equation for the 0(3)/ 0(2) nonlinear sigma 
model, whose solution is 

2 87T 

a = o"o = A exp — — (32) 

wherein A is a momentum cutoff. The fields an and <Jii must therefore be shifted by the 
constant o"o in order to fluctuate about zero. Then <7o multiplies |e 2 and in the large- 
N effective action; thus, it has the interpretation of a quantum fluctuation-induced mass 
squared for the fields e a . The excitations of these fields are then triplets with energy gap 
A = A exp — (4n/3g). After restoring the characteristic speed c, and the dependence on £, 
this becomes 



A = (Aa) J,i±i Vd - ^ - 2C + 1) exp - {f ^ fl( ^f - - l} } (33) 

Of course, this result neglects 1/N corrections. As £ — > 1 (the collinear point) the gap goes 
to zero. This signals a continuous phase transition, beyond which a different description of 
the system is needed, for instance in terms of the 0(3)/ 0(2) sigma model. From £ = 1 up to 
£ = 0, which corresponds to J' — ► oo, the gap increases and becomes eventually proportional 
to J'. 

The main feature of the excitations of this model is the degeneracy of the triplets. This 
ultimately comes from the SO (2) symmetry existing between the vectors ei and e2, which 
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may be mixed with each other without affecting the action. It is then not surprising if this 
degeneracy survives the large- N approximation. 

Indeed, 1/N corrections to the dispersion of the excitations can in principle be calculated, 
like they are in the 0(3)/ 0(2) sigma model (see, for instance, Polyakov's textS). The 
quantity of interest is then the self-energy of the fields e a . In general, this is a matrix 
T> a b(uj,k), which may be calculated diagrammatically with the help of the vertices shown 
in Fig. H] and of the propagators for the auxiliary fields a a b and <fi. The latter scale as 1/N 
since they are inverse polarization operators, and the order in 1/N is simply determined by 
the number of auxiliary field internal lines. However, the propagators for an and ct 22 are 
identical, since the large- ./V propagators of ei and e 2 are the same. Moreover, the structure 
of the vertices of Fig. |2| makes it impossible for a non-diagonal self-energy £ 12 to arise: 
even though the auxiliary fields ayi and may change a ei quanta into a e 2 quanta, the 
number of such vertices must be even, without tadpoles. The net result is that no self-energy 
diagram exists in which an entering ei line is turned into an exiting e 2 line, and consequently 
£12 = 0. The 1/N contributions to the En are show diagrammatically on Fig. Q At all 
orders in 1/N the two self-energies En and S 22 are identical and therefore the degeneracy 
of the excitations branches seems a exact feature of the model. 



IV. DISCUSSION 

In Sec II we argued that the frustrated AF spin chain may be described, in the continuum 
limit, by the 5*0(3) nonlinear sigma model defined by the Lagrangian fll7l) . In Sec III, we 
argued that the main feature of the spectrum of this continuum model is a doubly degenerate 
triplet of spin-1 excitations, separated from the singlet ground-state by a gap. However, this 
conclusion is incompatible with what is known of the spectrum of the frustrated spin-^ chain. 
The Lieb-Sc^Mattis (LSM) the0rem B 8tat es th- a sp in -| e hab h as e ither no g ap 0r 
degenerate ground states (corresponding to spontaneously broken parity). More recently, it 
has been shownEl that this theorem applies to half-integer spin chains in general. In the 
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special case J' = J/2 the exact ground state of the spin-| chain obtained by Majumdar and 
Gosh! is two-fold degenerate. This is compatible with the LSM theorem, but in contradiction 
with the field theory. Thus, the latter does not correctly describe frustrated half-integer spin 
chains. 

Before throwing the SO (3) model away, one should see if it describes the long wavelength 
behavior of integer spin chains. Exact results for the spin-1 chain are lacking and thus we are 
reduced to numerical study of finite systems. Quantum Monte-Carlo studies of frustrated 
spin systems are quite difficult because of a sign problem, therefore we limited ourselves 
to numerical diagonalizations of small chains (up to 14 sites). We used various methods, 
of which the most practical turned out to be a variant of the Lanczos methodS applied to 
subspaces of fixed momentum. Diagonalizations of spin-^ chains with up to 22 sites have 
also been performed, for the sake of comparison. It turns out that the spectra of spin- 
1 and spin-i chains of small lengths are qualitatively different. The ground state of the 
spin-1 chains is non-degenerate, with a gap to a doubly degenerate triplet of spin-1 states, 
exactly as in the field theory. Excitation spectra for J' = J/2 and various chain lengths 
are shown in Fig. |4], whereas spectra for 14 spins and various values of J'/ J are shown 
in Fig. |5|. The dispersion relation is symmetric with respect to parity (k — > —k) and the 
minimum of the spin-1 branch occurs at a wavenumber fc smaller than ir (in units of 1/a). 
Consequently, this minimum also occurs at — k , which is distinct from k except at the 
ferromagnetic (k = 0) and antiferromagnetic (k = 7r) points. The existence of these minima 
away from the parity-invariant points should not surprise us. If it were not for the quantum 
fluctuation disordering the system, there would be massless Goldstone branches around the 
ordering wave- vectors k = a and —k (a is given by Eq. (JD). The effect of quantum 
fluctuations is to produce a singlet ground state and to raise the minima of the Goldstone 
branches. We expect that the positions of these minima be approximately the same as the 
classical ordering wave- vectors ko and — fco, although we have no way to prove that they 
coincide exactly with these values. The 5*0(3) field theory effectively describes low-energy 
excitations about the minima of these excitation branches, and the degeneracy is essentially 
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linked to the fact that the classical ordering wavevectors ko and — ko are inequivalent. This 
feature of the excitation spectrum supports the view that the SO (3) field theory describes 
the long wavelength behavior of the frustrated spin-1 chain. 

It should be pointed out that, in a chain of length N, there is critical value of the 
ratio J/4J', above which the minimum in the triplet dispersion relation occurs at k = it. 
This critical value depends on the length of the chain. If one trusts the heuristic mapping 
of section II, it should be equal to ^ oc, - ) = 1 when N — > oo. For finite chains it is closer to 
|, but it is conceivable that it reaches the predicted value as N grows. Figure || shows the 
evolution of £ c for small sizes, as well as the value of the mass gap A. 

The question remains as to why this semi-classical treatment fails in the half-integer spin 
case. In other words, where in the argument of Sec. II should the distinction between integer 
and half-integer spin come into play? The answer may lie in our neglect of configurations with 
discontinuities. This has already been pointed out in Ref. (| in relation with the kinetic term. 
An easier way to see this is to consider the limit J' > J. In that limit one may reconsider 
the problem and view the system as a pair of interwoven, antiferromagnetic chains with a 
small interchain interaction J. Each chain could be semi-classically described by the usual 
0(3)/ 0(2) nonlinear sigma model, with the Lagrangian density (p, plus a topological term 
Ctop expressed as 

g 

£to P = ~2 e ' X dx< ^ ' ( 34 ) 

Since the effective lattice spacing for each chain is now a = 2a, the characteristic velocity 
would be c = 4J'as, while the coupling constant is still g = 2/s. One checks that this agrees 
perfectly with the £ — > limit of Eq. fl22|). Thus, as £ = J/4 J' decreases, the two unit 
vectors e x and e 2 of the SO (3) model become the fields of two weakly coupled 0(3)/ 0(2) 
sigma models, with the correct values of c and g. 

The interaction Ci nt between the two chains may be easily expressed in terms of ei and 
e 2 . Indeed, recall that in taking the continuum limit of a single AF chain, one uses the 
decomposition Sj = s[(— l) l e + 51], where 1 is the local magnetization and a = 2a is the 
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effective lattice spacing for each chain. Neglecting derivatives, the interaction Lagrangian 
may be written as 

L int = -2Js 2 ~a\ 1 -\ 2 . (35) 
The elimination of li and I2 is done by substituting the equation of motion0 

la = - jjr^( e a x d t e a ) . (36) 

4r O S(X 

The interaction Lagrangian then becomes 

Lint = ~Y^~J7^( ei X 5 * ei ) ' ( e2 X 9te2 ) ( 37 ) 

= — {ei ■ e 2 d t e 1 ■ d t e 2 - (e t ■ d t e 2 )(e 2 ■ d t ei)} . (38) 

If we assume that the two vectors ei and e 2 are identically orthogonal, as the classical 
ground state suggests if £ is small, then the above interaction exactly agrees with Eq. (|2T|) 
and the £ — >• limit of 7 as given in Eq. (p^). Moreover, the topological terms S , top [ei] and 
5'top[e2] cancel each other if ei is identically orthogonal to e 2 . This is particularly easy to 
understand if the configuration ei belongs to the homotopy class of the identity, since the 
orthogonality constraint then effectively makes e 2 a mapping from the sphere 5*2 to the circle 
Si and 7r 2 (Si) = 0. In general, the two topological terms cancel simply because of the fact 
that two orthogonal unit vectors specify a rotation and 7r 2 (SO(3)) = 0. Thus, if we assume 
that ei _L e 2 in the small £ limit, the Lagrangian (|2l|) is recovered with exactly the same 
parameters, since the topological terms are absent. 

However, the orthogonality ei _L e 2 is not strict, but only favored energetically. A local 
deviation from this orthogonality gives back their full importance to the topological terms 
and causes a distinction between integer and half-integer spins. Such a distinction does not 
occur in the 5*0(3) formulation of the problem since the orthogonality is then 'built-in'. In 
that formulation, a local deviation from orthogonality corresponds to a discontinuity in the 
order parameter, which was not allowed from the start. 
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Finally, let us point out that the above remarks concerning parity-breaking ordering 
wavevectors in frustrated antiferromagnets and their consequences on the excitation spec- 
trum may apply to other systems, such as the nearest-neighbor Heisenberg antiferromagnet 
on a triangular lattice. In this case it is believed that the ground state has long-range order. 
However, it is conceivable that the introduction of vacancies destroys this order without 
affecting the applicability of the SO (3) theory: the effect of the vacancies would then be 
to increase the value of the effective coupling constant g beyond the critical value g c . A 
priori, it is not clear if the SO (3) field theory describes the long wavelength behavior of 
antiferromagnets on the Kagome lattice or on the see-saw chain0. A semi-classical study of 
these systems may be interesting. 
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FIGURES 




FIG. 1. Classical ground state of the frustrated AF chain with J'/ J = \. Each spin makes a 
120° angle with its neighbors and the periodicity is the smallest possible, i.e., three sites. Below is 
different look at the same spin chain, using a 'railroad trestle' geometry. 




FIG. 2. Vertices used in the l/N expansion of the model (27). 




FIG. 3. l/N contributions to the self-energy En of the field ei. 
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FIG. 4. Excitation spectrum of the frustrated spin-1 chain with J' = J/2 for chains lengths 
L = 8, 10, 12 and 14. The white circles represent spin-triplet states, while the black squares 
represent spin-singlets, when lower in energy. 
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FIG. 5. Excitation spectrum of the 14-site frustrated spin-1 chain with J' / J = 1, \ and \. 
Again, the white circles represent spin-triplet states and the black squares spin-singlets, when lower 
in energy. 
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FIG. 6. Dependence on chain length iV of the critical coupling ratio £ c = (J/4J') C at which the 
triplet degeneracy appears (error bars). Also illustrated is the energy gap A/J for ^ = ^ (circles). 
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